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ABSTRACT 

Symbiotic binaries are comprised of nebulae, whose densest portions have electron 
concentrations of 10 8 — 10 12 cm~ 3 and extend to a few AU. They are optically thick 
enough to cause a measurable effect of the scattering of photons on free electrons. 
In this paper we introduce modelling the extended wings of strong emission lines by 
the electron scattering with the aim to determine the electron optical depth, r e , and 
temperature, T c , of symbiotic nebulae. We applied our profile-fitting analysis to the 
broad wings of the Ovi 1032, 1038 A doublet and Hen 1640 A emission line, measured 
in the spectra of symbiotic stars AG Dra, Z And and V1016 Cyg. Synthetic profiles fit 
well the observed wings. By this way we determined r e and T c of the layer of electrons, 
throughout which the line photons are transferred. During quiescent phases, the mean 
t c = 0.056 ± 0.006 and T c = 19 200 ± 2 300K, while during active phases, mean 
quantities of both parameters increased to r c — 0.64 ± 0.11 and T e = 32 300 ± 2 000K. 
During quiescent phases, the faint electron-scattering wings are caused mainly by free 
electrons from/around the accretion disk and the ionized wind from the hot star with 
the total column density, N c < 10 23 cm~ 2 . During active phases, the large values of r 
are caused by a supplement of free electrons into the binary environment as a result 
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> . 1 INTRODUCTION 



Symbiotic stars are long-period interacting binaries with 
orbital periods in the range of years. In their spectra we 
can recognize three main sources of radiation. The first 
one is represented by a giant star of spectral type (G- 
)K-M , the second one is a very hot (21 > 10 5 K) com- 
pact star, which is in most cases a white dwarf accret- 
ing from the wind of the giant. The third component 
of radiation is produced by a nebula, which represents 
the ionized fraction of the c i rcumstellar material in the 
binary (e.g. Boyarchuk 1967; Scaquist, Taylor & Button 



ing effects of Rayleigh, Raman and Thomson scattering. Ra- 
man and Rayleigh scattering results from interaction be- 
tween the hot star photons and the neutral atoms in the gi- 
ant's wind. They are important tools in mapping the ioniza- 
tion structure of symbiot i c binaries (e.g 
iNussbaumer et al. Ill989l; ISchmid \\l998 : 



Isliker et al"lll989l; 



Birriel 



2004; [Lee I 



1984 hereafter STBl; iKenyon I [l986| ; INussbaumer fc Vogel 



1984, ncrcaltcr SIB; Kcnyon iy 
19871 ; ICorradi et al. 1 120031 ; ISkop 



aal 1 120051 . and references 
therein). As a result, the circumstellar environment of sym- 
biotic binaries comprises energetic phot ons from the hot star 
with luminosities of 1 2 - 10* (e.g. iMiirset et al~lll99ll ; 
iGreiner et al. I IT997I ; ISkopal I l2005t ). neutral particles pro- 
duced by th e cool giant at rates of a few times 10 ~ 7 Mp yr _1 
(e.g. STB, iMikolaiewska et al. I l2002i ; ISkopal 1 120051 ). and 
ions and free electrons resulting from the processes of ioniza- 
tion. Symbiotic stars thus represent ideal objects for study- 
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20091 . |2012h ." In contrast, the Thomson scattering of photons 
by free electrons acts within the ionized part of the symbiotic 
stars environment, and thus can diagnose the symbiotic neb- 
ula. In the spectrum, we can indicate this effect in the form 
of shallow, wide wings of the strongest emission lines, which 
photons are scattered by free electrons, and thus are Doppler 
shifted by their thermal motion to both the red and blue side 
of the line. The effect of this process is weak and wavelength 
independent, because of a very small and constant value of 
the Thomson cross-section, <jt = 6.652 x 10~ 25 cm 2 . Never- 
theless, the densest portions of symbiotic nebulae with elec- 
tron concentrations of log(n c ) ~ 8 — 12 (n e in cm -3 ), could 
be optically thick enough to cause a measurable effect of the 
electron scattering. From this point of view, strong emission 
lines of highly ionized elements that are formed in the dens- 
est part of the ionized medium in a vicinity of the hot white 
dwarf, represent the best candidates. 

Originally, ISchmid et all (l999) suggested that the 
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broad wings of the Ovi 1032 and 1038 A resonance lines 
could be expl ained by scattering of the Ovi photons by 
free electrons. lYoung et al. 1 |2005) successfully compared a 
model of the electron scattering wings with the Ovi dou- 
bl et in the AG Dra s pectrum. This process was also used 
by I Jung fc Leel l|2004) to model the broad Ha wings in the 
spectrum of the symbiotic star V1016 Cyg as an alternative 
to t he Raman scatterin g of Ly/3 photons on atomic hydro- 
gen. ISkopal et alj l|2009f ) used the electron scattering of the 
Ovi 1032, 1038 A doublet and Hen 1640 A line in the AG Dra 
spectrum to support the origin of the X-ray-UV flux anti- 
correlation, revealed by modelling the SED. 

In this paper we describe a simplified method for fitting 
the broad wings of intense emission lines by the electron 
scattering process (Sect. 3). We applied our profile-fitting 
procedure to the Ovi resonance doublet and the Hen 1640 A 
line, observed in the spectra of symbiotic stars AG Dra, 
Z And and V1016 Cyg, with the aim to determine the elec- 
tron optical depth, r e , and temperature, T c , of their nebulae. 
Sect. 4 presents the results and their discussion. Conclusions 
are found in Sects. 5. 



2 OBSERVATIONS AND DATA TREATMENT 

For the purpose of this paper, we used extremely in- 
tense emission lines of the Ovi 1032, 1038 A doublet ob- 
served in the FUSE (Far Ultraviolet Spectroscopic Explorer), 
BEFS (Berkeley Extreme and Far-UV Spectrometer) and 
TUES (Tubingen Ultraviolet Echelle Spectrograph) spectra 
of AG Dra, Z And and V1016 Cyg. We used also the strong 
emission line Hen 1640 A in the high-resolution IUE (Inter- 
national Ultraviolet Explorer) spectra of AG Dra. The spec- 
tra were obtained from the satellite archives with the aid 
of the Multimission Archive at the Space Telescope Science 
Institute (MAST). They are summarized in Table 1. 

The FUSE spectra were processed by the calibration 
pipeline version 3.0.7, 3.0.8 and 2.4.1. We used the calibrated 
time-tag observations (TTAG photon collecting mode). Be- 
fore adding the flux from all exposures we applied an ap- 
propriate wavelength shift relative to one so to get the best 
overlapping of the absorption features. Then we co-added 
spectra of individual exposures and weighted them accord- 
ing to their exposure times. The wavelength scale of the 
spectra was calibra ted with the aid of the i nterstellar ab- 
sorption lines (e.g. iRogerson fc Ewell 1 1 1984] ). Accuracy of 
such calibration is of ± 0.05 A. Finally, we binned the re- 
sulting spectrum within 0.025 A. 

All the spectra were corrected for heliocentric ve- 
locity including that of the satellite, and dered dened 

" Il99ll ). 0.28 
for AG 



Table 1. Log of spectroscopic observations. 



with E B -v = 0-30 (Z And. iMurset et al 
(V10 16 Cyg. iNussbaumer fc Schild Ill98ll 1 and 



Dra dBirriel et al. I 120001 1, 
ICardelli et al. I j 19891 ). 



using the extinction curve of 



3 THE MODEL 

3.1 Assumptions and simplifications 

Thomson scattering represents a special case of the scat- 
tering of a photon off a free electron, which is fully elas- 
tic (i.e. the photon energy does not change). In the real 



Date 


Julian date 


Orb. 


Obs. 


Line'') 


(dd/mm/yyyy) 


- 2 400 000 


phase" ) 


Satellite 




AG Dra 


11/12/1981 


44949.5 


0.38 


IUE 


Hen 


18/09/1993 


49248.5 


0.22 


BEFS 


Ovi 


29/06/1994 


49532.5 


0.74 


IUE 


Hen 


09/07/1994 


49542.5 


0.75 


IUE 


Hen 


12/07/1994 


49545.5 


0.76 


IUE 


Hen 


28/07/1994 


49561.5 


0.79 


IUE 


Hen 


17/09/1994 


49612.5 


0.88 


IUE 


Hen 


28/07/1995 


49927.5 


0.46 


IUE 


Hen 


14/02/1996 


50127.5 


0.82 


IUE 


Hen 


22/11/1996 


50409.5 


0.33 


TUES 


Ovi 


16/03/2000 


51619.5 


0.54 


FUSE 


Ovi 


25/04/2001 


52024.5 


0.28 


FUSE 


Ovi 


14/11/2003 


52957.5 


0.98 


FUSE 


Ovi 


24/06/2004 


53180.5 


0.38 


FUSE 


Ovi 


25/12/2004 


53364.5 


0.72 


FUSE 


Ovi 


15/03/2007 


54174.5 


0.20 


FUSE 


Ovi 


Z And 


05/07/2002 


52460.5 


0.90 


FUSE 


Ovi 


04/08/2003 


52855.5 


0.42 


FUSE 


Ovi 


V1016 Cyg 


10/08/2000 


51766.5 


? 


FUSE 


Ovi 



a ) - According to iFekel et alH feOOd) , 
1032, 1038 A doublet, 



b ) - Hen 1640 A, Ovi 



case, the photon always transfers some of its energy to the 
electron, which shifts its wavelength by ~+0.024 A, the so- 
called Compton wavelength of the electron. However, the 
Doppler effect arising from the thermal motion of free elec- 
trons leads to significantly larger shifts to both sides of 
the spectrum. Therefore, the elastic Thomson scattering is 
a good approximation in studying scattering of low-energy 
photons (hv « m e c 2 ) off non-relativistic electrons (see e.g. 
iRosswog fc Briiggen 1120071 . in detail). 

To model the effect of the electron scattering we 
have to know how the scattered photons are redis- 
tributed in frequenc i es an d directions. It was shown by 
iHummer fc Mihalas I (| 1967t > that it is convenient and suf- 
ficiently accurate to regard the radiation field, from which 
scattering occurs, as isotropic, so that the direction can be 
averaged out. Therefore, for the sake of simplicity, we con- 
sider isotropic scattering with Maxwellian distribution of 
electron velocities. Under these assumpti ons the redistri bu- 
tion function can be expressed in a form dMihalas 1119701 ). 



R, 



(v, V ) = — 



e 

V / 



2w 



erfc 



2w 



(1) 



where v and v are frequencies of radiation before and after 
the scattering and w is the electron Doppler width, 



c 



2kT 



(2) 
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and the complementary error function erfc(x) is defined as 

2 



erfc (x) 



dz. 



(3) 



To calculate the electron-sc attering wings profiles, we 
adopted a simplified scheme of I Munch I (|l950l ) , which as- 
sumes that a plane-parallel layer of free electrons of the op- 
tical thickness r c and the temperature T c is irradiated by 
the line photons, and that the electrons are segregated from 
the other opacity sources, which implies no change in the 
equivalent width of the line. In our case, this assumption 
corresponds to modelling the wings of highly ionized Ovi 
and Hen lines, which are formed within the +5 and/or 
He + zone close to the hot white dwarf photosphere, and the 
layer of free electrons, where the Thomson scattering arises, 
is located above the line formation region. 

According to the radiative transfer equation, assuming 
that the electron scattering is the only process attenuating 
the original line photons, the observed line flux can be ex- 
pressed as 



F° 



Foe" 



(4) 



where Fo is the line flux before scattering, r c = otA^ is the 
electron optical depth and N c is the column density of free 
electrons along the line of sight. As the scattered fraction 
of the original flux is redistributed into the line wings so 
that the equivalent width before and after the scattering is 



constant, F w 



= Fo - F obs = F (l 



Then the line 



profile after it emerges from the la yer of scattering elec trons 
may be approximated by (see also lCastor et al. II1970T ) 



To $(a;) + (1 - e~ Tc ) / <f>(x )R e (x ,x)dx , (5) 



where is the incident line profile, R c (x , x) is the redis- 
tribution function for Thomson scattering (Eq. (1)) and x or 
x is a frequency displacement from the line center in units 
of electron Doppler width before and after the scattering, re- 
spectively. So, the first right-side term of Eq. represents 
the original flux at x attenuated by the scattering (Eq. @), 
and its scattered, (1 — e _T °), fraction is redistributed in the 
wings (the second term) . Note that the scattered profile used 
previously by e.g. ICastor et al. I (|l970l ) represents a special 
case of Eq. © for r e -C 1. 

To fit the theoretical profile (fSJ) to observations, means 
to determine its variables, r c , T c , and those of the incident 
profile $(k). 



3.2 The incident profile 

First, we estimated the continuum level from a large wave- 
length interval around the Ovi lines by a linear fit to the 
noise in the spectrum. In the case of the Hen 1640 A line, 
the continuum level was estimated with the aid of the corre- 
sponding low-resolution spectrum. Second, we approximated 
the incident profile $(A) in the model (J3J) as follows. The Hen 
1640 A line was possible to fit with a single Gauss curve. The 
fit of the observed emission core provided the first estimate 
of its position, width and the height. The Ovi 1032 A line 
is in most cases asymmetric with respect to its reference 
wavelength. Its blue emission wing is steeper than the red 
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Figure 1. Examples of the observed (solid lines) and the incident 
(dotted lines) profiles of the Ovi 1032 A line. Arrows denote the 
red-side emission bumps. 



one, being cut by an absorption component (see Fig. [!}. It 
is probably caused by the scattering in the line at the close 
vicinity of the white dwarf within the densest part of the 
wind moving to the observer. This absorption can be more 
pronounced duri ng active pha ses, when a higher mass-loss 
rate is observed (|Skopalll200"r3 ). The observation of AG Dra 
from 2007 May 15, made during the 2006-08 active phase, 
is consistent with this view (see the left panel of Fig. [T]). 
The scattering in the line operates within the line forma- 
tion region. Therefore, we take the incident profile of the 
Ovi 1032 A line as a sum of emission and absorption Gaus- 
sians. The incident emission of the Ovi 1038 A line was not 
possible to reconstruct by a direct fitting of its observed re- 
mainder, because of a strong influence by th e absorption of 
the in terstellar molecular hydrogen, H2 (e.g. ISchmid et al. I 
1999). Therefore, we reconstructed the incident Ovi 1038 A 
profile with the aid of the theoretical ratio of the doublet 
lines with the assumption that they have the same width 
(see below). 

Based on the above mentioned observational properties 
of the Ovi 1032, 1038 A doublet, we reconstructed its inci- 
dent profile by a superposition of three Gaussians, 



$(A) = ^/ n exp 



A - X r 



(6) 



where indices n = 1 and 3 denote the Gaussians of the emis- 
sion cores at Ai ~ 1032 A and A 3 ~ 1038 A. The second 
curve (n = 2) represents the absorption component, which 
cut the blue side of the 1032 A emission line (Fig. [1} . How- 
ever, this rather strong absorption did not reproduce the 
narrow one seen in the scale of Figs. [2] and [3] as a P-Cyg 
component in the 1032 A line, but not recognizable in the 
scale of Fig. [1] Only for the V1016 Cyg spectrum, due to a 
symmetrical emission core of the 1032 Aline, it was possible 
to fit this narrow P-Cyg absorption by the I2 component 
(see the right panel of Fig. [3}. Irrespectively of its origin, we 
neglected its influece to the profile. 

The effect of the interstellar H2 absorption to the 
Ovi 1038 A line profile was significant in the spectra of Z And 
and V1016 Cyg, because of a large amount of interstellar 
matter on the line of sight to these objects (Fb-v ~ 0.3). 
As a result, the observed intensity ratio I1/I3 was w 7 for 
Z And and m 15 for V1016 Cyg. Therefore, we reconstructed 
the original Ovi 1038 A line adopting the theoretical ratio 
I1/I3 = 2, <T3 = o\ and A3 so to fit the emission remainder 
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of the line. The absorption effect from the H2 molecules was 
apparently fainter on the AG Dra spectra (Eb-v ~ 0.08, 
Ii /I3 ~ 2) , and thus the estimate of the initial Ovi 1038 A 
line profile was more trustworthy. 

Some spectra show a noticeable additional emission fea- 
tures at the red side of the 1032 A and 1038 A line cores up 
to ~ 10% of their height. We did not investigate the ori- 
gin of these features and did not take them into account 
in the fitting procedure. However, in the AG Dra spec- 
trum from 2007 March, the unknown emission features were 
rather strong. Therefore, we included them to the incident 
Ovi 1032 A, 1038 A line profiles as an additional Gaussians 
(see Fig. [T|. 



T e reflect different properties of the symbiotic nebula dur- 
ing different levels of the activity. During quiescent phases, 
the mean T c = 19 200 ± 2 300K, r c = 0.056 ± 0.006, while 
during active phases T c = 32 300±2000K, r c = 0.64±0.11, 
respectively. Here the uncertainties represent rms errors of 
the average values. 

The average quantities of T e and r e as derived from 
quiescent and active phases agree well wit h those found 
independently by modelling the SED (e.g. ISkopal 1 120051 ; 
ISkopal et al. 2009). More than one order of magnitude dif- 
ference in t c (i.e. in N c ) reflect a significant increase of free 
electrons on the line of sight in the direction to the hot star 
during active phases. The origin of this change is shortly 
discussed in Sects. 4.2. and 4.3. 



3.3 Profile-fitting analysis 

First, we selected the flux-points of the observed profile(s), 
F obs (Ai), for fitting with the function ((5]). They were selected 
from the observed profiles by omitting some artificial (sharp) 
emission/absorption features and the depression around the 
Ly/3 line, caused by the Rayleigh scattering (see Fig. [2}. To 
find the best solution, we calculated a grid of models for rea- 
sonable ranges of the fitting parameters, 7 n , a n , A n for the 
original profile and r c , T e for the scattered wings. In 

the case of the Ovi 1032 A, 1038 A doublet, the parameters 
73, (73, A3 of its 1038 A component were estimated accord- 
ing to the properties of both the doublet lines, as described 
in Sect. 3.2. The grid of models was prepared with steps 
Ar c = 0.0005, AT C = 100 K, A7 n = 7„/100, A<r„ = 0.001 
and AA n = 0.001 A. Finally, we selected the model corre- 
sponding to a minimum of the function 



2 

Xred 



1 N 

-T 



F obs (Ai) - *(A ; ) 



AF obs (Ai) 



(7) 



where F° bs (Ai) are the observed fluxes of the profile, TV is 
their number (» 200), N is the number of d.o.f., AF obs (Ai) 
are their errors and ^(Ai) are theoretical fluxes. 

Errors in the selected flux points were around of 10- 
15% of the line wings. Based on them we determined un- 
certainties in T c and r c for individual spectra. To obtain 
a rough estimate of the corresponding range of T c we 
fixed other fitting parameters and varied T e so to fit fluxes 
F° bs (Ai)±AF obs (Ai). Similarly we proceeded to estimate the 
limits for r c . Such the estimated ranges of fitting parameters 
can be as large as ~ 50% of the best model value and are 
asymmetrically placed with respect to it (see Table 2). This 
is a result of the non-uniformly distributed fluxes determin- 
ing the profile of electron-scattered wings (see the begin- 
ning of Sect. 3.3.). In the case of the AG Dra spectrum from 
14/02/1996, the relatively small uncertainties in T c and r e , 
but the large value of \rcd probably reflect underestimated 
flux errors, AF obs (Ai). Large flux uncertainties on the spec- 
trum from 22/11/1996 and poorly defined continuum level 
did not allow us to estimate the upper limit of T (see Fig. [2]). 



4 RESULTS AND DISCUSSION 

The resulting fits are shown in Fig. [2] and [3] and correspond- 
ing parameters are in Table 2. Significant changes in r e and 



4.1 Application to selected symbiotic stars 

4.I.I AG Dra 

AG Dra is a yellow sy mbiotic star comprising a cool giant of 
a K2 III spectral type jMurset fc Schmid ll999T ). It is located 
at a high galactic latitude of 41°, which implies that its 
spectrum is less affected by the interstellar matter. As a 
result, the line ratio 7(1038 A)/7(1032 A) was close to its 
theoretical value of 0.5, which made the modelling of the 

Ovi doublet more trustworthy (Sect. 3.2). 

By modelling the UV/IR continuum of AG Dra. lSkopal I 

(2005) found that the mean electron temperature of the neb- 
ula during quiescent phase and/or small bursts runs between 
18 000 and 21 800 K, while during major outbursts (1980-81, 
1994-5, 2006-7, Fig. [2]), the nebula significantly strengthens 
and increases its mean T e to ~ 35 000 K. In this work we con- 
firmed these results independently by the profile-fitting anal- 
ysis of the electron scattered wings. We found that during a 
quiescent phase, the mean T c = 21 700 ± 3 600 K, while dur- 
ing active phases T e = 32 300 ± 2 000 K. Also r c is a function 
of the star's activity. Our analysis revealed r c = 0.063±0.007 
and 0.64 ± 0.11 during quiescence and activity, respectively 
(Table 2, Fig. [2] and . 

In the 15/03/2007 spectrum, a relatively high value of 
r c = 0.20, but very faint wings of the Ovi doublet (see 
Fig. [2J is caused by the weak incident line flux, Fo, be- 
cause F w i n g/Fo = 1 - e~ T ° (Sect. 3.1., Fig. [Q. So, the well 
detectable electron scattering wings reflect a relatively large 
amount of free electrons on the line of sight during the tran- 
sition from the major 2006 outburst. 



4.1.2 Z And 

Z And is considered to be a prototype of symbiotic stars. 
Here the white dwarf accretes from the stellar win d of a 
M3-4III red giant fe.g. lFernandez-Castro et al. Ifililsl ). Dur- 
ing active phases, the light curve shows 2-3 mag bright- 
cnings, while the quiescent phase is ch aracterized by a 
wave-like orbitally-related variation (e.g. Belvakina~l 1 19851 ; 
iFormiggini fc Leibowitz Ill994l ; ISkopal et al. 1120061 ). 

Two FUSE spectra used in this work were observed 
at the end of the major 2000-03 outbursts, 05/07/2002, 
and just after the optical re brightening on 04/08/2003 (see 
Fig. 1 of lSkopal et al. ||2006| ). The electron scattering wings 
observed during both dates (Fig. [3} corresponded to very 
similar quantities of the fitting parameters, T a ~ 16 500 K, 
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Figure 2. Top p anel shows the U ligh t curve of AG Dra from 1977. Active phases are characterized by outbursts with multiple maxima. 
The data are from lSkopal et aL Bottom panels compare the observed (dotted + enhanced line) and modelled (solid smooth line) 

broad wings of the Ovi 1032, 1038 A doublet and the Hen 1640 A line at different stages of activity. The enhanced parts of the observed 
profile were fitted with the function J5j. Horizontal dotted line represents the level of the continuum. Timing of individual observations 
are given in Table 1 (arrows in the top panel). Fluxes are in units of 10 — 12 erg cm -2 s _1 A -1 . 
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Figure 3. As in Fig. [2] but for Z And and V1016 Cyg. 



Table 2. Best solutions of Eq. 10 for r c and T a and their ranges Ar c and AT e , corresponding to AF obs (Ai) (see Sect. 3.3). 



Date 


Stage* 






T c 


AT 




(dd/mm/yyyy) 








[K] 


[K] 




AG Dra 


11/12/1981 


A 


0.39 


0.26-0.55 


40100 


22 500-47 200 


0.6 


18/09/1993 


T 


0.18 


0.15-0.21 


20 200 


14 600-23 100 


1.9 


29/06/1994 


A 


0.39 


0.34-0.44 


32 700 


29 700-34 900 


1.6 


09/07/1994 


A 


1.17 


1.0-1.4 


31000 


28 800-34 500 


3.6 


12/07/1994 


A 


0.86 


0.70-1.1 


29100 


25 300-38 700 


3.0 


28/07/1994 


A 


0.80 


0.69-0.95 


41300 


40 100-46 700 


2.5 


17/09/1994 


A 


0.49 


0.36-0.66 


25 400 


14 900-27 400 


1.2 


28/07/1995 


A 


0.39 


0.33-0.46 


28 800 


27500-33 500 


0.9 


14/02/1996 


T 


0.30 


0.26-0.35 


14 300 


13 700-15 800 


12.1 


22/11/1996 


T 


0.23 


0.13-0.34 


23 900 


18 600->90 000 


11.6 


16/03/2000 


Q 


0.044 


0.033-0.056 


12 800 


7600-18 500 


1.0 


25/04/2001 


Q 


0.065 


0.051-0.078 


23 100 


22 900-25 200 


1.9 


14/11/2003 


T 


0.084 


0.073-0.096 


29 900 


16 200-41 900 


1.2 


24/06/2004 


Q 


0.063 


0.051-0.074 


20 300 


11200-24 000 


0.8 


25/12/2004 


Q 


0.078 


0.066-0.091 


30 700 


29 100-38 900 


1.8 


15/03/2007 


T 


0.20 


0.16-0.25 


21300 


9 600-27 900 


2.3 


Z And 


05/07/2002 


Q 


0.050 


0.042-0.058 


16 500 


8 600-22 300 


2.3 


04/08/2003 


Q 


0.051 


0.045-0.056 


15 300 


10 300-22 400 


1.5 


V1016 Cyg 


10/08/2000 


Q 


0.038 


0.031-0.044 


16 000 


8 200-26 100 


2.0 



* A - active phase, Q - quiescent phase, T - transition to quiescence 



t c ~ 0.050 and T c ~ 15 300 K, r c ~ 0.051, respectively (Ta- 
ble 2). Electron temperatures are consistent with t hose de- 
rived from modelling the SED in the continuum ()Skopal I 
l2005h . 

4.1.3 V 1016 Cyg 

V1016 Cyg is a member of a small group of symbiotic stars 
called symbiotic novae. In 1964 it underwent a nova-like out- 
burst jMcCuskev II 19651 ), during which the star's brightness 
increased from m pg ~ 15.5 to ~ 10.5 mag in 1971, following 
a gradual small decrea se to ~ 11.7 in 2000 (see Fig. 1 of 
iParimucha et al. 1120021 ). 

There was only one well exposed spectrum in the FUSE 
archive, from 10/08/2000. The wings of the Ovi doublet are 
clearly visible, although the contamination of the 1038A line 
is significant (Eb-v = 0.28 mag). In this spectrum we were 



able to fit also the sharp absorption component in the P- 
Cyg profile of the 1032 A line (Fig. [3|. Fitting parameters 
of the nebula, r c = 0.038 and T c = 16 000 K are similar to 
those derived from observations during quiescent phase of 
AG Dra. The very small value of r c is a result of a very 
weak wings with respect to the strong central emission core, 
F wing /Fo = 0.037. 

4.2 Thomson scattering during quiescent phase 

Fitting the extended faint wings of the Ovi and Hen lines 
by Thomson scattering suggested a connection between r e 
and the level of the star's activity (Fig. [2j. 

During quiescent phase, it is assumed that the sym- 
biotic nebula arises from ionizing a portion of the neutral 
wind from the giant only, i.e. the wind from the hot star 
is neglected. As the source of neutral particles and that of 
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ionizing photons are separated, the nebula will be spread 
asymmetrically around the hot star in the binary. This im- 
plies that the column densities of free electrons on the line 
of sight to the hot star, and thus r e , will also be a func- 
tion of the orbital phase, (p. Here we investigate the function 
r e (ifi) for the simplest (idealized) case as outlined by STB. In 
particular, we assume a spherically-symmetric unperturbed 
wind from the giant, whose particles are accelerated along 
the /3-law, and the stationary situation (i.e. no binary rota- 
tion and no gravitational attraction of the accretor to the 
wind are included). Under these assumptions, the extent of 
the ionized zone during quiescence can be obtained from a 
parametric equation 



X = f(r,4), 



(8) 



which solution defines the boundary between neutral and 
ionized gas at the orbital plane, determined by a sys- 
tem of polar coordinates, r, with the origin at the hot 
star. The function /(r, #) was treated for the first time by 
STB for a steady sta t e situ ation and pure hydrogen gas. 
INussbaumer fc Vogel I (| 19871 ) considered also a contribution 
of free electrons from singly ionized helium, because its zone 
nearly overlaps that of the Hn region. This increases the 
electron concentration in the ionized zone by a factor of 
(l + o(He)), where a(He) is the abun dance by number of He 
relativ e to H. Following derivation of INussbaumer fc Voeel I 
( 1987), but replacing the terminal velocity of the wind, «oo, 
by its /3-law distribution 



Uwind(r) =v ao (l- — ~\ 



(9) 



where r is counted from the centre of the cool giant with the 
radius R s (= the beginning of the wind), we can express the 
terms of Eq. @ as 



X 



Q B (H,T c )(l + a(He)) 



pLh 



Ah 



(10) 



and 



0.5 - 



- 



-0.5 - 



Giant 



Hot 
star 



D 



HI 



HII 



-0.5 



0.5 
Distance / p 



1.5 




Figure 4. Top panel shows the ionization structure of AG Dra 
as seen pole-on. The Hl/Hn boundary was calculated according 
to Eq. JSJ and X = 8 (Sect. 4.2.). Bottom panel displays the 
corresponding function r c (</?) for the orbital inclination of 60 and 
90°. Labels A,B,C,D mark directions at specific orbital phases. 
Crosses correspond to r c = 0.044—0.078, derived from the spectra 
during quiescent phase (Table 2). 



2/3 



du. 



(11) 



The ac celeration param eter in the wind ([9]), /3 = 2.5 for red 
giants (|Schroder I fl985). /i is the mean molecular weight, 
77ih is the mass of the hydrogen atom, qb(H, T e ) stands for 
the total hydrogenic recombination coefficient for case B, p 
is the separation of the binary components, Lu is the flux 
of hydrogen ionizing photons (s~ ) and M g is the mass-loss 
rate from the giant. Finally, we expressed the radial dis- 
tances s and r in units of p, defining u = s/p and I = r/p = 
Vu 2 + 1 - 2u cos 0. Solutions of Eq. ((8| for ug(= sg/p) at 
directions 6 define the ionization Hi/Hn boundary. For the 
AG Dra parameters (p = 355 Rq, Lu = 3.4 x 10 46 s _1 , 



lo- 



33 R^-,, Mr, 



13 cm 3 s" 1 



3.2 x IP' 7 Mf^yr' 1 
~aLTll995l : 



Mikolaicwska et 



Woo = 30 kms 
and qb = 1-4 x 

iFekel et al. 1 120001 : ISkopaT 2005; Skopal et al. 200i) the ion- 
ization parameter X = 8, which corresponds to an open Hn 
zone around the hot star (see Fig. 3| . 

Having defined the ionization structure in the binary, we 
can calculate the function r c (ip) by integrating the electron 
concentration throughout the ionized zone along the line of 



sight to the hot star, i.e. 



T e = ITT 



3 (s)ds 



(12) 



where the electron concentration n e (s) = (1 + o(He))n(r) 
and n(r) is the density of hydrogen atoms in the wind from 
the giant. If the line of sight is not passing through the 
neutral region, sg — > 00 (in practice we adopted sg — 100 x 
p). To calculate correctly T c (i p) for AG Dra, we cons idered 
its orbital inclination i = 60° (|Schmid fc Schild Ill997h . The 
function t c (i/p) is plotted in Fig. [4] A maximum value of t c 
(~0.03) is around ip — 0.1, when the line of sight passes 
the ionization region as an asymptote to the boundary (the 
direction B in the figure). A minimum of ~0.007 corresponds 
to the position with the hot star in front (ip — 0.5), because 
of the lowest densities of the wind from the giant. Figure [4] 
demonstrates that theoretical T c (cp) function is significantly 
below the values derived from observations during quiescent 
phases (crosses in the figure, Table 2). This implies that the 
ionized fraction of the unperturbed wind from the giant is 
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not capable of producing the observed electron scattering 
wings. 

In a more realistic case the density distribution in a 
binary with mass-losing giant is determined mainly by the 
rotation of the binary and accretion by its compact com- 
panion, as w as demostrated by several hydrodynamical sim- 
ulations (e.g. Theuns fc Jorissen ^993; Bisi kalo et al. |[l995l ; 
iMastrodemos fc Morris Ill998l ; iNagae et al. 112004 ). In these 
studies, it was shown that the regions with highly increased 
density around the accretor, around the mass-losing star and 
the whole binary, and behind the accretor (opposite to its or- 
bital motion) were in the form of a disc, a spiralling arm and 
an elongated accretion wake, respectively. As a result, the 
column density of free electrons in the direction to the ac- 
cretor can be enriched by the ionized material accumulated 
at/around the accretion disk at each position of the binary. 
Additional extremes can be expected for highly inclined 
orbits, when the line of sight pass es throughout a hig her 
density structure (see e.g. Fig. 6 of iDumm et al. ifeoOOl ). A 
lower orbital inclination smooth s out the density contrast s 
at different phases (see Fig. 2 of lTheuns fc Jorissen |[T993h . 
Accordingly, in the real case, the faint electron-scattering 
wings during quiescent phases with r c = 0.044 — 0.078 can 
be caused mainly by free electrons from/around the accre- 
tion disk and the ionized wind from the hot star, which 
corresponds to the total column density, N c = t c /gt < 
10 23 cm -2 . The presence of the latter was proved by more 
autho rs (e.g. IVogel Ill993r . iNussbaumer et al. |[l995l ; ISkopal I 

l20oeh . 



4.3 Thomson scattering during active phase 

It is well known that during active phases the hot stars 
in symbi otic binaries enhance significantly the mass-loss 
rate ( e.g. iFernandez-Castro et ai~ 1995 : INussbaumer et al. I 



ll995l ; ICrocker et al. Iboo'l " ISkopal 1120061 ). The ejected mate- 
rial is ionized by the luminous central hot star, which thus 
enh ances ra d iation from the symbiotic nebula. For exam- 
ple, ISkopal I ((2005) derived a factor of «10 stronger nebu- 
lar emission in the continuum during active phases with re- 
spect to values measured during quiescence. The enhanced 
mass-loss rate from the hot star thus represents a significant 
supplement of free electrons into the nebula. As a result, 
the electron optical depth, r c = <7t-/V c , will be considerably 
larger during active phases (Table 2). Here, this is well doc- 
umented by the series of the AG Dra spectra (Fig. The 
model SED demonstrates that the nebular radiation domi- 
nates also the spectral region of the photometric U filter (see 
top panels of Fig. [SJ. Therefore, the level of the AG Dra ac- 
tivity is well mapped with the U light curve (top panel of 
Fig- E]), which thus explains the relationship between r e and 
the star's brightness in U (bottom panel of Fig. [SJ . A large 
scatter around r c = 0.2 is probably connected with the tran- 
sition phase, when the nebula can be partially optically thick 
in the continuum. 

We can conclude that during active phases the large 
values of r c = 0.39 - 1.17 (i.e. N c = 0.58 - 1.8 x 10 24 cm" 2 ) 
are caused by a supplement of free electrons into the binary 
environment as a result of the enhanced mass-loss rate from 
the hot star. 




3.2 3.3 3.4 3.5 

log(Wavelength) [A] 



E 
3 













O 








$ o 








o 










o 


o 


o 










o 










3 0.1 


0.2 


0.3 


0.4 0.5 


0.6 0.7 0.8 0.9 1 


1.1 1.2 



Figure 5. Top panels show the IUE spectra of AG Dra during 
the 1994 active phase (SWP51632+LWP28752) and a quiescent 
phase (SWP37473+LWP16675). The dashed, solid thin and dot- 
ted line represent the continuum from the nebula, the hot star and 
the giant, resp ectively. The solid thick line is their sum, the model 
SED (see ISkopal ll2005l . in detail). The gray line figures schemat- 
ically the profile of the photometric U filter. Bottom panel shows 
the dependence of r e on the [/-magnitude (Fig. [2] Table 2). 



5 CONCLUSIONS 

We investigated the effect of Thomson scattering of the 
strong emission lines, Ovil032A, 1038 A and Hen 1640 A, 
observed in the spectra of symbiotic stars AG Dra, Z And 
and V1016 Cyg. Our models of their profiles are in good 
agreement with those observed by FUSE, BEFS, TUES and 
IUE satellites (Figs. [2] and [3]). This supports the idea that 
the broad wings of these lines result from the scattering of 
the line photons on free electrons. Their profile is given by 
two fitting parameters, r c and T c (Table 2), which charac- 
terize the scattering region, i.e. the symbiotic nebula. Par- 
ticular results of our analysis can be summarized as follows. 

(i) The presence of the electron scattering wings in the 
line profiles of highly ionized elements locates their origin 
mostly to the vicinity of the hot star in the binary, with the 
highest density on the line of sight. 

(ii) During quiescent phases, the mean T e = 19 200 ± 
2 300K, while during active phases T c = 32 300 ± 2 000K. 
This findings agrees well with quantities derived indepen- 
dently by modelling the SED. 

(iii) The electron optical depth also depends strongly on 
the star's activity. During quiescent phases, r c = 0.056 ± 
0.006, while during active phases, r c = 0.64 ±0.11 (Table 2, 



Parameters of symbiotic nebulae from Thomson scattering 9 



Figs. [2]and[3]|. Uncertainties in the results given in (ii) and 
(iii) represent rms errors of the average values. 

(iv) During quiescent phases, the ionized fraction of the 
wind from the giant within a simple (STB) model (Sect. 4.2) 
is not capable of giving rise a measurable effect of the Thom- 
son scattering. In the real case, the faint electron-scattering 
wings with r e = 0.044 — 0.078 are caused mainly by free elec- 
trons from/around the accretion disk and the ionized wind 
from the hot star with the total column density of a few 
xlO 22 cm -2 (Sect. 4.2., Fig.©. 

(v) During active phases, the large values of r c are caused 
by a supplement of free electrons into the nebula from the 
enhanced wind of the hot star, which increases N c to ~ 
10 24 cm" 2 (Sect. 4.3., Pig-EJ. 

The presumable relationship between r c and the level of the 
activity, as suggested by our profile-fitting analysis, could 
be used in probing the mass-loss rate from the hot stars in 
symbiotic binaries. 
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